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Abstract 


This paper considers the problem of controlling a class of 
nonlinear multibody flexible space systems consisting of a flexible 
central body to which a number of articulated appendages are 
attached. Collocated actuators and sensors are assumed, and 
global asymptotic stability of such systems is established under 
a nonlinear dissipative control law. The stability is shown to 
be robust to unmodeled dynamics and parametric uncertainties. 

For a special case in which the attitude motion of the central 
body is small, the system, although still nonlinear, is shown to 
be stabilized by linear dissipative control laws. Two types of linear 
controllers are considered: static dissipative (constant gain) and 
dynamic dissipative. The static dissipative control law is also 
shown to provide robust stability in the presence of certain classes 
of actuator and sensor nonlinearities and actuator dynamics. The 
results obtained for this special case can also be readily applied for 
controlling single-body linear flexible space structures. For this 
case, a synthesis technique for the design of a suboptimal dynamic 
dissipative controller is also presented. The results obtained in 
this paper are applicable to a broad class of multibody and single- 
body systems such as flexible multilink manipulators, multipayload 
space platforms, and space antennas. The stability proofs use 
the Lyapunov approach and exploit the inherent passivity of such 

systems. 

1. Introduction 

A class of spacecraft envisioned for the future will require flexible multibody systems such 
as space platforms with multiple articulated payloads and space-based manipulators used for 
satellite assembly and servicing. Examples of such systems include the Earth Observing System 
(EOS) and Upper Atmospheric Research Satellite (UARS). (See refs. 1 and 2, respectively.) 
These systems are expected to have significant flexibility in their structural members as well as 
joints. Control system design for such systems is a difficult problem because of the hi ghl y 
nonlinear dynamics, large number of significant elastic modes with low inherent damping, 
and uncertainties in the mathematical model. The published literature contains a number of 
important stability results for some subclasses of this problem (e.g. , linear flexible structures, 
nonlinear multibody rigid structures, and most recently, multibody flexible structures). Under 
certain conditions, the input-output maps for such systems can be shown to be passive. (See 
ref. 3.) The Lyapunov and passivity approaches are used in reference 4 to demonstrate global 
asymptotic stability (g.a.s.) of linear flexible space structures (with no articulated appendages) 
for a class of dissipative compensators. The stability properties were shown to be robust to first- 
order actuator dynamics and certain actuator-sensor nonlinearities. Multibody rigid structures 
represent another class of systems for which stability results have been advanced. Subject 
to a few restrictions, these systems can be ideally categorized as natural systems. (See ref. 5.) 
Such systems are known to exhibit global asymptotic stability under proportional-and-derivative 
(PD) control. After recognition that rigid manipulators belong to the class of natural systems, 
a number of researchers (refs. 6—9) established global asymptotic stability of terrestrial rigid 
manipulators by using PD control with gravity compensation. Stability of tracking controllers 
was investigated in references 10 and 11 for rigid manipulators. In reference 12, the results 
of reference 11 were extended to accomplish exponentially stable tracking control of flexible 



multilink manipulators local to the desired trajectory. Lyapunov stability of multilink flexible 
systems was addressed in reference 13. However, the global asymptotic stability of nonlinear, 
multilink, flexible space structures, which to date has not been addressed in the literature, is 
the principal subject of this paper. 

The structure to be considered is a complete nonlinear rotational dynamic model of a 
multibody flexible spacecraft which is assumed to have a branched geometry, (i.e., a central 
flexible body with various flexible appendage bodies). (See fig. 1.) Throughout this paper, 
actuators and sensors are assumed to be collocated. Global asymptotic stability of such systems 
controlled by a nonlinear dissipative controller is proved. In many applications, the central 
body has large mass and moments of inertia when compared with that of the appendage bodies. 
As a result, the motion of the central body is small and can be considered to be in the linear 
range. For this case, robust stability is proved with linear static as well as with linear dynamic 
dissipative compensators. The effects of realistic nonlinearities in the actuators and sensors are 
also investigated. The stability proofs use Lyapunov’s method (ref. 9) and LaSalle’s theorem. For 
systems with linear collocated actuators and sensors, the stability proof by Lyapunov’s method 
can take a simpler form if the work-energy rate principle (ref. 13) is used. However, because the 
work-energy rate principle is applicable only when the system is holonomic and scleronomic in 
nature, a more direct approach is used here in evaluating the time derivative of the Lyapunov 
function which makes the results more general. 

The organization of the paper is as follows. In section 3, a nonlinear mathematical model of 
a generic flexible multibody system is given. A special property of the system, which is pivoted 
to the proofs, and some basic kinematic relations of the quaternion (i.e., measure of attitude of 
the central body) are given. Section 4 establishes the global asymptotic stability of the complete 
nonlinear system under a nonlinear control law based on quaternion feedback. A special case, in 
which the central body attitude motion is small, is considered in section 5. Global asymptotic 
stability is proved under static dissipative compensation, and these results are extended to 
the case in which certain actuator and sensor nonlinearities are present. In addition, dynamic 
dissipative compensators, a more versatile class of dissipative compensators, are considered. 
Section 6 presents the application of the results to the important special case of linear single- 
body spacecraft. Two numerical examples are given in section 7 and some concluding remarks 
are given in section 8. 

2. Symbols 


A 

system matrix 

A 

diagonal matrix (eq. (79)) 

B 

control influence matrix 

B 

diagonal matrix (eq. (80)) 

C 

Coriolis and centrifugal force coefficient matrix 

C 

diagonal matrix of squares of eigenvalues 

D 

damping matrix 

D 

structural damping matrix 

E 

Euler transformation matrix; modulus of elasticity 

F 

compensator output matrix 
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plant transfer function 
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actuator transfer function 
acceleration gain matrix 
position gain matrices 
rate gain matrix 
estimator gain matrix 
k x k identity matrix 
index variable 

inertia matrix of single-body spacecraft 

performance index 

index variable 

stiffness matrix of system 

stiffness matrix of system for linear case 

stiffness matrix for flexible degrees of freedom 

controller transfer function 

number of rigid body degrees of freedom 

Lagrangian of system 

extended Lebesgue space 

mass-inertia matrix of system 

mass-inertia matrix of system for linear case 

number of total degrees of freedom 

number of flexible degrees of freedom 

Lyapunov matrix 

generalized coordinate vectors 

vector of rigid body coordinates 

output weighting matrices 

vector of flexural coordinates 

control weighting matrix 

skew-symmetric matrix 

property of the system 

Laplace variable 

terminal time; kinetic energy 

transfer matrix 

vector of control input 

Lyapunov function candidate 

integrator output 


3 



state vectors 
acceleration output 
position output 
rate output 
state vector 
the quaternion 

ith component of unit vector along eigenaxis 
vector part of quaternion 
ith component of quaternion 
scalar defined as ( a 4 — 1 ) 
integral of uj 
transformation matrix 
performance function 
state vector 
Euler angle vector 

vector of rotational degrees of freedom between rigid bodies 
diagonal matrix (eq. (82)) 
scalar gain 
input to nonlinearity 
state vector 

damping in ith structural mode 
controller state vector 
mode shape matrix 
Euler rotation 
defined as —9 

actuator nonlinearity (ith loop) 
position sensor nonlinearity (ith loop) 
rate sensor nonlinearity (ith loop) 
cross product operator of vector u> 
angular velocity vector for central body 
Abbreviations: 

a.s. asymptotically stable 

DDC dynamic dissipative controller 

EOS Earth Observing System 

g.a.s. globally asymptotically stable 
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LQG linear quadratic Gaussian 

MSPR marginally strictly positive-real 

PD proportional and derivative 

PR positive-real 

SPR strictly positive-real 

UARS Upper Atmospheric Research Satellite 

3. Mathematical Model 

3.1. Equations of Motion 

The class of systems considered herein consists of a branched configuration of flexible bodies 
as shown in figure 1. Each branch by itself could be a serial chain of structures. For the sake of 
simplicity and without loss of generality, consider a spacecraft with only one such branch where 
each appendage body has one degree of freedom (hinge joint) with respect to the previous body 
in the chain. The results presented in this paper are also applicable to the general case with 
multiple branches. Consider the spacecraft consisting of a central flexible body and a chain of 
( k — 3) flexible links. The central body has three rigid rotational degrees of freedom, and each 
link is connected by one rotational degree of freedom to the neighboring link. The Lagrangian 
for this system can be shown (ref. 14) to have the following form: 

L = p T M{p)p-q T Kq (1) 

where p = 9^ , q^)^ , w is the 3x1 inertial angular velocity vector (in body-fixed coordi- 

nates) for the central body, 9 = (#i, 02, . . . , 8(k-3)) T (^i denotes the joint angle for the zth joint 
expressed in body-fixed coordinates), q is the (n - k ) vector of flexible degrees of freedom (modal 
amplitudes), M{p) = M^{p) > 0 is the configuration-dependent mass-inertia matrix, and K is 
the symmetric positive-definite stiffness matrix related to the flexible degrees of freedom. From 
the Lagrangian equation (1), the following equations of motion axe obtained: 

M(p)p + C(p,p)p + Dp + Kp = B T u (2) 

where p = (l T ,0 T ,q T ) T , 7 = w, C(p,p) corresponds to Coriolis and centrifugal forces, D is the 
symmetric positive-semidefinite damping matrix, B = [hxk °jfcx(n-jfc) ] is the control influence 
matrix, and u is the k vector of applied torques. The first three components of u represent the 
attitude control torques applied to the central body about its X- , Y -, and Z-axes; the remaining 
components are the torques applied at the ( k - 3) joints. The symmetric positive-semidefinite 
stiffness and damping matrices K and D are 

Ofcxfc ^kx(n-k) 1 f Ofcxfc Ofcxfn— k) 

K= _ D= ; (3) 

-®(n—k)xk K(n—k)x(n—k). .^(rt-k)xk ^(n-k)x(n-k). 

where K and D are symmetric positive-definite. The details of the derivation of the mathematical 
model can be found in reference 14; a summary is provided in the appendix. 

The angular measurements for the central body are Euler angles (not the vector 7), whereas 
the remaining angular measurements between bodies are relative angles. One important inherent 
property of such systems that is crucial to the stability results presented in this paper is defined 
next. 
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Property S: The matrix (\M - C) is skew-symmetric for the system represented by 
equation (2). The justification of this property is given by theorem A1 in the appendix. 

The central-body attitude (Euler angle) vector 77 is given by £( 77)17 = w, where £( 77 ) is a 3 x 3 
transformation matrix. (See ref. 15.) The sensor outputs consist of three central-body Euler 
angles, the ( k — 3) joint angles, and the angular rates (i.e., the sensors are collocated with the 
torque actuators). The sensor outputs are then given by 

y P = Bp y T = Bp (4) 

where p = {p T 1 0 T , q T ) T and 77 is the Euler angle vector for the central body. The measured 
angular position and rate vectors are y p = (77 ,9 T ) T and y r = (u> T ,9^)^, respectively. The 
body rate measurements ui are assumed to be available from rate gyros. The input-output map 
from u to y r is passive as shown in theorem A2 of the appendix. 


3.2. Quaternion as Measure of Attitude 

The orientation of a free-floating body can be minimally represented by a three-dimensional 
orientation vector. However, this representation is not unique. Euler angles are commonly used 
as a minim al representation of the attitude. As stated previously, the 3x1 Euler angle vector 77 
is given by £( 77)77 = oj, where £( 77 ) is a 3 x 3 transformation matrix. For certain values of 77 , 
£( 77 ) becomes singular; however, note that the limitations imposed on the allowable orientations 
because of this singularity are purely mathematical in nature and have no physical significance. 
The problem of singularity in the three-parameter representation of attitude has been studied 
in detail in the literature. An effective way of overcoming the singularity problem is to use the 
quaternion formulation. (See refs. 16-19.) 

The unit quaternion (Euler parameter vector) a is defined as follows: 



(5) 


where 5 = (Si, S 2 , &z) T is the unit vector along the eigenaxis of rotation and <p is the magnitude 
of rotation. The quaternion is subject to the norm constraint 

a T a + a\ = 1 ( 6 ) 


The quaternion can also be shown (ref. 19) to obey the following kinematic differential equations: 


o 


-(cj x a + 0:401) 


(7) 


04 = 


1 T— 
- -or o 

2 


( 8 ) 


The attitude control of a single-body rigid spacecraft with quaternion feedback has been 
thoroughly investigated in references 16-19. The quaternion representation is used here for 
the central-body attitude. The quaternion can be computed from Euler angle measurements 
given by equations (4). (See ref. 20.) 

The open-loop system, given by equations (2), (7), and (8), has multiple equilibrium solutions 
(oJ s ,O 4 5S , 0 ^) T , where the subscript ss denotes the steady-state value; the steady-state value 
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of q is zero. By defining /? = (q 4 - 1) and denoting p = 2 , equations ^2), (7), and (8) can be 
rewritten as v ' 

M z + C z 4- Dz + Kq — B^u 
= [ 0 (n-3)x3 -^(n— 3)x(n— 3) ] z 

O = X a +(/? + 1)0;] 

In equation (9a) the matrices M and C are functions of p and (p,p), respectively. Note that the 
farst three elements of p associated with the orientation of the central body can be fully described 
by the unit quaternion. Hence, M and C are implicit functions of a, and therefore the system 

represented by equations (9)-(ll) is time-invariant and can be expressed in the state-space form 
as follows: 

x = f(x,u) (J 2 ) 

where x = (a T , 0, 6 T , q T , z T ) T . Note that the dimension of x is (2 n + 1), which is one more 
t an the dimension of the system in equation (2). However, the constraint of equation (6) is 
now present. Verification that the constraint of equation (6) is satisfied for all t > 0 if it is 
satisfied at t = 0 easily follows from equations (7) and (8). 

4. Nonlinear Dissipative Control Law 

Consider the dissipative control law u, given by 


9 

? 


(9a) 

(9b) 

( 10 ) 

( 11 ) 


U - -G p p - G r y r ( 13 ) 

where p = (a T , 0 T ) T . Matrices G p and G r are symmetric positive-definite k x k matrices- G n is 
given by ’ p 

(l + (^— 03x(fc-3) 

° r= „ r (14) 

L U (A-3)x3 G P2 (fc _ 3)x{A ._ 3) _ 

Note that equations (13) and (14) represent a nonlinear control law. If G p and G r satisfy certain 
conditions, this control law can be shown to render the time rate of change of the system’s energy 
negative along all trajectories (i.e., it is a dissipative control law). 

The closed-loop equilibrium solution can be obtained by equating all the derivatives to zero 
in equations (2), (10), and (11). In particular, p = p = 0 => u> = 0,9 = 0,q = 0, and 


-B t G p p 


-G p p 


0/fcxl 

,°(n-fc)xl. 


.Kq . 


(15) 


Because ofequation (6), \(3 + 1| < 1. Therefore G p is positive-definite and equation (15) implies 
p = {a ,6 ) T = 0 and q = 0. The equilibrium solution of equation (11) is 0 = 0 SS = Constant 
(i.e., a 4 = Constant), which implies from equation (6) that a 4 = ±1. Thus, there appear to 
be two closed-loop equilibrium points corresponding to a 4 = 1 and a 4 = -1 (all other state 
variables are zero). However, from equations (5), a 4 = 1 ^ = 0, and a 4 = -1 => <f> = 2n (i.e., 

only one equilibrium point exists in the physical space). 

One of the objectives of the control law is to transfer the state of the system from one orien- 
tation (i.e., equilibrium position) to another orientation. Without loss of generality, the target 
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orientation can be defined to be zero, and the initial orientation, given by"(a(0), 04(0), 0(0)), 
can always be defined in such a way that |<fo(0)| < 7r, 0 < 04(0) < 1 (corresponds to \<j>\ < x), 
and (a(0),c*4(0)) satisfy equation (6). 

The following theorem establishes the global asymptotic stability of the physical equilibrium 
state of the system. 

Theorem 1: Suppose G p2 {k _ 3]x() ._ 3] and G rfcxfc are symmetric and positive-definite, and 
G p i = nh, where p > 0. Then, the closed-loop system given by equations (12) and (13) is 
globally asymptotically stable (g.a.s.). 

Proof: Consider the candidate Lyapunov function 

V = \i> T M<j>)p + l -q T Kq 4- ±9 T G p2 9 + \a T {G pl + 2pl 3 )a + p0 2 (16) 

Here, V is clearly positive-definite^ and radially unbounded with respect to a state vector 
(a T ,/?,0 T , q T ,p T ) T because M(p), K, G p 1, and G p2 are all positive-definite symmetric matrices. 
Note that the matrix M(p), although configuration dependent, is uniformly bounded from below 
and above by the values which correspond to the minimum and maximum inertia configurations, 
respectively (i.e., there exist positive-definite matrices M and M such that M < M < M). The 
time derivative of V results in 

V = p T Mp -1- l -p T Mp + q T Kq + 9 T G p2 9 + ^(Gpi + 2 pl 3 )a -I- 200pt (17) 

With the use of equations (2), (4), (10), (11), and (14), 

V = p T B T u + - C^jp - p T Dp - p T Kp + q T Kq + 9 T G p2 9 

+ i(fia) r G p ia + i(/3 + l)u; r Gpia + pu T a (18) 

where Q. = (wx) denotes the skew-symmetric cross product matrix (i.e., uxi = fix). With the 
substitution for u, the fact that p T Kp = q T kq and (fta) T G pl a = 0, and the use of property S 
of the system, equation (18) becomes 

V — —p^^D + GrB^p — {Bp) T G p p + — (/? 1 )w'^Gpio -1- puj^ oc + 9^ G p2 9 (19) 

Note that ( Bp) T G p p = %(0 4- l)w T G p ia + pu T a + 9 T G p2 9. After several cancellations, 

V = -p T (d + B T G r B')p (20) 

Because ( D + B T G r B) is a positive-definite symmetric matrix, V < 0 (i.e., V is negative- 
semidefinite) and V r = 0=^p = 0=>p = 0. By substitution in the closed-loop equation, equa- 
tion (15) results. As shown previously, equation (15) => p = 0 and q = 0; i.e., a = 0, 0 = 0, and 
<*4 = ±1 (or 0 = 0 or — 2). Consistent with the previous discussion, these values correspond to 
two equilibrium points representing the same physical equilibrium state. 

Prom equation (16), verification easily follows that any small perturbation e in 04 from the 
eq uili brium point corresponding to 04 = —1 will cause a decrease in the value of V (e > 0 because 
< 1). Thus, in the mathematical sense, 04 = —1 corresponds to an isolated equilibrium 
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point such that V = 0 at that point and V < 0 in the neighborhood of that point (i.e., 04 = — 1 
is a repeller and not an attractor). Previously, V has been shown to be negative along all 
trajectories in the state space except at the two equilibrium points. That is, if the system’s 
initial condition lies anywhere in the state space except at the equilibrium point corresponding 
to 0:4 = —1, then the system will asymptotically approach the origin (i.e., x = 0); if the system 
is at the equilibrium point corresponding to 04 = -1 at t = 0, then it will stay there for all t > 0. 
However, this is the same equilibrium point in the physical space; hence, by LaSalle’s invariance 
theorem, the system is globally asymptotically stable. Q.E.D. 

5. Systems in Attitude-Hold Configuration 

Consider a special case where the central-body attitude motion is small. This can occur i 
many real situations. For example, in cases of space station-based or shuttle-based manipulators, 
the inertia of the base (central body) is much greater than that of any manipulator link 
or payload. In such cases, the -rotational motion of the base can be assumed to be in the 
linear region, although the payloads (or links) attached to it can undergo large rotational 
and translational motions and nonlinear dynamic loading because of Coriolis and centripetal 
accelerations. The attitude of the central body is simply 7 (the integral of the inertial angular 
velocity ui) and the use of quaternions is not necessary. The equations of motion (2) can now be 
expressed in the state-space form simply as 

- f ° 1 

X ~[-M~ 1 K + D) 

where x = {jp T ,p T ) T , p = (■y T , 0 T , q T ) T , and M and C are functions of x. 

5.1* Stability With Static Dissipative Controllers 

The static dissipative control law u is given by 

u — Gpyp G r yr (22) 

where G p is symmetric positive-definite k x k matrix and 

Vp = Bp y r = Bp (23) 

where y p and y r are measured angular position and rate vectors, respectively. 

Theorem 2: Suppose c Pfcxt and Gr kxk are symmetric and positive-definite. Then, the 
closed-loop system given by equations (21)-(23) is globally asymptotically stable. 

Proof: Consider the candidate Lyapunov function 

V = i p T M(p)p + (k + B T G p B^jp (24) 

Clearly, V is positive-definite because M(p) and (K + B T G p B ) are positive-definite symmetric 
matrices. Defining K = (K + B T G p B), the time derivative of V can be shown to be 

V = p T Qm - c'jp -p T Kp +p T Kp-p T (d + B T G T B^p (25) 

Again, with the use of property S, p T {\M - C)p = 0, and after some cancellations, 

V= -p T (p + B T GrB)p (26) 


0 

b t 


( 21 ) 
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Because (D + B T G r B) is a positive-definite matrix, V < 0 (i.e., V is negative-semidefinite in p 
and p and U = 0=>p = 0=>p = 0). With substitution in the closed-loop equation (2) and u 
given by equation (22), 

(K + B T G p B^p = 0 => p = 0 (27) 

Thus, V is not zero along any trajectories; then, by LaSalle’s theorem, the system is g.a.s. Q.E.D. 

The significance of the two results presented in theorems 1 and 2 is that any nonlinear 
mnl t.i hod v system belonging to these classes can be robustly stabilized with the dissipative 
control laws given. In the case of manipulators, this means that any terminal angular position 
can be achieved from any initial position with guaranteed asymptotic stability. 


5.2. Robustness to Actuator-Sensor Nonlinearities 

Theorem 2, which assumes linear actuators and sensors, proves global asymptotic stability for 
systems in the attitude-hold configuration. In practice, however, the actuators and sensors have 
nonlinearities. The following theorem extends the results of reference 4 (pp. 59—62) to the case 
of nonlinear flexible multibody systems. That is, the robust stability property of the dissipative 
controller is proved to hold in the presence of a broad class of actuator-sensor nonlinearities 
with the following definition: a function ip{v) is said to belong to the (0, oo) sector (fig. 2(a)) if 
ip(0) = 0 and vrp{v) > 0 for v ± 0 and ip is said to belong to the [0, oo) sector if vip{v) > 0. 

Let i p a i{-), tjjpi(-), and Wt(') denote the nonlinearities in the zth actuator, position sensor, 
and rate sensor, respectively. Both G p and G r are assumed to be diagonal with elements Gpi 
and G r i, respectively; then the actual input is given by 

Ui — tpai \~Gjriippi (ypi) ~~ G r i‘tp r i{y r i^ (i = 1,2, k) (28) 

With the assumption that ip pi, tpai, and ip ri (z = 1,2, .. . ,k) are continuous single-valued func- 
tions, R — *■ R, the following theorem gives sufficient conditions for stability. 

Theorem 3: Consider the closed-loop system given by equations (21), (23), and (28), 
where G p and G r are diagonal with positive entries. Suppose xp a i, ip pi, and ip r i are single-valued, 
time-invariant continuous functions, and that (for z = 1,2, . . . ,k) 

1. ip a i are monotonically nondecreasing and belong to the (0,oo) sector 

2. xppi and ip ri belong to the (0, oo) sector 

Under these conditions, the closed-loop system is globally asymptotically stable. 

Proof: (The proof closely follows that of ref. 4.) Let ip = -y p {k vector). Define 

4>pi ( 29 ) 

‘ t Pri(v)= ( 3 °) 

liippi, ipri e (o, oo) or [0, oo) sector, then z/^, tp ri also belong to the same sector. Now, consider 
the following Lure-Postnikov Lyapunov function: 


k y 

V = \p T M{p)p + \q T Kq + Y] f tpai [^pi ^ pt ‘( I/ )] dv 
2 2 i=I Jo 


(31) 
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where K is the symmetric positive-definite part of K. Differentiation with respect to t and use 
of equation (2) yield 


k 

V = P t (b t u - Cp-Dp - Kp) + '-p T Mp + £ 1 + fkg 

i = 1 

Upon several cancellations and the use of property S, 

k k 
V = E UlVri ~ ViH’ai 

i=l i=l 

where matrix D is the positive-definite part of D. 



(32) 


(33) 


(34) 


Because are monotomcally nondecreasing and tp ri belong to the (0, oo) sector, V < 0 and 
the system is at least Lyapunov-stable. In fact, it will be proved next that the system is g a s 
inrst, consider a special case when ijj ai are monotonically increasing. Then, V < -g T Dg and 

V ~ 0 only w ? en * = 0 ^d <P = 0, which implies p = 0 => p = 0. Substitution in the closed-loop 
equation results in 

Kp = B T xp a [-GptPp (j/p)] (35) 

0 _ rfa [ — Gpippfyp}] 

Kg J ~ [ 0 J (36) 

^ 'tpa [~ — 0 and ^ = 0 

Ifippi belong to the (0, oo) sector, tp a i( u ) = ippiiv) = 0 only when v = 0. Therefore y p = 0. Thus. 
^ = 0 only at the origin, and the system is g.a.s. 

In the case when actuator nonlinearities are of the monotonically nondecreasing type (e.g. 
saturation nonlinearity), V can be 0 even if Cp # 0. Figure 2(b) shows a monotonically 
nondecreasing nonlinearity. However, every system trajectory along which V = 0 will be shown 
to go to the origin asymptotically. When <p # 0, V = 0 only when all actuators are locally 
saturated. Then, from the equations of motion, the system trajectories will go unbounded, 
w ich is not possible because the system was already proved to be Lyapunov-stable. Hence, 

cannot be identically zero along the system trajectories, and the system is g.a.s. Q.E.D. 

For the case when the central-body motion is not in the Unear range, the results of robust 
stability in the presence of actuator-sensor nonlinearities cannot be easily extended because the 
stabiUzing control law given in equation (13) is nonlinear. 

The next section extends the robust stabihty results of section 5.1 to a class of more versatile 
controllers caUed dynamic dissipative controUers. The advantages of using dynamic dissipative 
controllers include higher performance, more design freedom, and better noise attenuation. 

5.3. Stability With Dynamic Dissipative Controllers 

To obtain better performance without the loss of guaranteed robustness to unmodeled 
dynamics and parameter uncertainties, consider a class of dynamic dissipative controllers (DDC) 
Such compensators have been suggested in the past for controlfing only the elastic motion 
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(refs. 21-23) of linear flexible space structures with no articulated appendages (i.e., single- 
body structures). These compensators were based on the fact that the plant, which consisted 
only of elastic modes with velocity measurements as the output, was passive (i.e., the transfer 
function was positive-real (PR)). A stability theorem based on Popov’s hyperstability concepts 
(ref. 24) was then used, which essentially states that a positive-real system controlled by a 
strictly positive-real (SPR) compensator is stable. Even in the linear single-body setting, certain 
problems occur with these results. In the first place, no attempt is made to control the rigid- 
body attitude, which is the main purpose of the control system. Secondly, the results assume 
that measurements of only the elastic motion are available and that the actuators affect only 
the elastic motion. These assumptions do not hold for real spacecraft unless the actuators are 
used in a balanced configuration for accomplishing only damping enhancement with no rigid- 
body control. (See ref. 4.) Finally, the stability theorem invoked assumes the compensator 
to be strictly positive-real, which is overly restrictive. It is also an ambiguous term having 
several nonequivalent definitions. (See ref. 25.) In view of these facts, the concept of marginally 
strictly positive-real (MSPR) systems will be introduced, which is stronger than the standard 
positive-real systems but is weaker than the weakly strictly positive-real systems defined in 
reference 25. 

The results of this section address the problem of controlling both rigid and elastic modes; 
these results essentially extend and generalize the results of reference 26, which also addressed the 
control of rigid-plus-flexible modes, but only for the linear single-body case. In this section, the 
stability of dynamic dissipative compensators for flexible nonlinear multibody space structures in 
the attitude-hold configuration is proved by using some of the results and methods of reference 26. 

5.3.1. Mathematical preliminaries. PR and MSPR systems are defined in definition 1 
and definition 2, respectively. 

Definition 1: A rational matrix- valued function T(s ) of the complex variable s is said to be 
positive-real if all of its elements are analytic in Re[s] > 0 and T(s) + T t (s*) > 0 in Re[s] > 0, 
where an asterisk (*) denotes the complex conjugate. 

Scalar PR functions have a relative degree (i.e., the difference between the degrees of the 
denominator and numerator polynomials) of —1, 0, or 1. (See ref. 27.) Positive-real matrices can 
be shown to have no transmission zeros or poles in the open right half of the complex plane, can 
have only simple poles on the imaginary axis with nonnegative definite residues. By application 
of the maximum modulus theorem, it is sufficient to check for positive-semidefiniteness of T{s ) 
only on the imaginary axis s = ju>, 0 < u < oo; i.e., the condition becomes T(ju) + T*(ju>) > 0. 
Suppose (A, B, C, D) is an nth-order minimal realization of T(s). From reference 28, a necessary 
and sufficient condition for T(s) to be positive-real is that there exists an n x n symmetric 
positive-definite matrix P and matrices W and L such that 


a t p + PA = -LL t ' 

c = b t p + w t l > 
w T w = d + d t . 


(37) 


This result is generally known in the literature as the Kalman-Yakubovich lemma. A stronger 
concept along these lines is the SPR systems. However, as stated previously, there are several 
definitions of SPR systems. (See ref. 25.) The concept of weakly SPR (ref. 25) appears to be 
the least restrictive definition of SPR. Nevertheless, all the definitions of SPR seem to require 
the system to have all poles in the open left half plane. 

The concept of marginally strictly PR systems is defined in reference 29 and included here 
in definition 2. 
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Definition 2: A rational matrix-valued function T(s) of the complex variable s is said 
to be marginally strictly positive-real (MSPR) if T(s) is PR and if T(ju) + T*(ju ) > 0 for 
u € (-00,00). 

The obvious difference between this definition and the definition of PR systems is that the 
weak inequality (>) has been replaced by strict inequality. The difference between the MSPR 
and weak SPR of reference 25 is that the latter definition requires the system to have poles in the 
open left half plane, whereas the former definition permits poles on the jui- axis. Reference 29 
shows that an MSPR controller can robustly stabilize a PR plant. 

5.3.2. Stability results. Consider the system given by equation (21) with the sensor outputs 
given by equations (23); then 

Vp ~ (7 r , T y T = (u T , 9 T ) T (38) 

where y p and y r axe the measured angular position and rate vectors, respectively. The 
central-body attitude rate measurements u> are assumed to be available from rate gyros. 

Suppose a controller £(s), with k inputs and k outputs, is represented by the minimal 
realization 


x c — A c x c + B c u c 


(39) 


Vc "b 


(40) 


where x c is the n c -dimensional state vector, (A c , B c , C c , D c ) is a minimal realization of Us) 
and u c = y p . 


Define 


Equations (39)-(43) can be combined as 


where 


A, = 


A c 0 
C c 0 


II 


(41) 

f T T'] 
z = (Xc . v ) 

r 

(42) 

II 


(43) 

— A z x z + B z u c 

(44) 

bt 

II 


(45) 

D c 

c,=[o 4] 

(46) 


The closed-loop system is shown in figure 3. The nonlinear plant is stabilized by fC(s ) if the 

closed-loop system (with fC(s) represented by its minimal realization) is globally asymptotically 
stable. 
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Theorem 4: Consider the nonlinear plant in equation (21) with y p as the output. Suppose 
that 

1. The matrix A c is strictly Hurwitz 

2. An (n c + k) x (n c + k) matrix P z = Pj > 0 exists such that 

A t z P z + P Z A Z = -Q z = -diag(ijLc,O fc ) (47) 

where L c is the k x n c matrix such that ( L c , A c ) is observable, and L c (sl — A c ) l B c has 
no transmission zeros in Re[s] > 0 

3. C z = BjP z (48) 

4. The controller IC(s) = C c (sl - A c )~ l B c + D c has no transmission zeros at the origin 
Then K(s) stabilizes the nonlinear plant. 

Proof: First consider the system shown in figure 4(a). The nonlinear plant is given by 
equation (2), and its state vector is taken to be {q*,p T ) T (i.e., {'l T ,0 T ) T is not included in the 
state vector). Now consider the Lyapunov function 

V = \p T M(j>)p + \q T Kq + \x T z P z x z (49) 

where K is the symmetric positive-definite part of K (i.e., the part associated with nonzero 
stiffness). Note that V is positive-definite in the state vector {q^ ,p T )^ because the mass-inertia 
matrix M(p), as stated previously, is symmetric positive-definite and uniformly bounded from 


below and above. Then 

V = p r M{p)p + ^ p T Mp + q T Kq + ^ ( x^P z x z + x T z P z x z ) (50) 

After substitution for M{p)p from equation (2) and for x z from equation (44), equation (50) 
becomes 

V = p T B T u — q T Dq + p T — c'jp — p T Kp + q T Kq 

+ | ^x t z A t z -t- ujB'[ > jP z x z + x t z P z {A z x z + B z u c )] (51) 

Property S of the system makes the matrix {\M — C ) skew-symmetric , and 

V = p T B T u - q T Dq + \x T z { A?P i + P x A z )x z + i £(bJp,)x, + \x T z (P z B z ) u e (52) 

V = -q T Dq + p T B T u — ^x^Q z x z + x^Cju c (53) 

where equations (47) and (48) have been used. Noting that u = — y z = —C z x z and Bp = t/r = w c 
(fig. 4(a)), 

V = -q T Dq - -x T z Q z x z - u*y z + yju c (54) 

from which t 

V = -q T Dq - -X T Z Q Z X Z (55) 
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Because D is positive-definite, V < 0 (i.e., V is negative-semidefinite) and the system is 
Lyapunov-stable. Now, V = 0 only if g = 0 and L c x c = 0. Therefore, either y r = 0 or y r 
consists only of terms such as ut e 2 o f where v is a constant vector and z 0 is a transmission 
zero of (A C ,B C ,L C ). No transmission zeros of ( A C ,B C ,L C ) exist in Re(s) > 0, which requires 
that y r —> 0 exponentially. Because (A c , B c , L c ) is minimal and stable, x c — ► 0 exponentially. 
But, y T — *■ 0 => 9 — > 0 and u> — > 0 =*• p — ► 0; this then implies that p — > 0. By substitution in 

equation (2), 7 7 ss,9 0 SS , q -> 0, and u-*0, where 7 SS and 6 SS are some steady-state 

values of 7 and 0, respectively. 

Now consider the configuration shown in figure 4(b), which is realized by application of the 
following similarity transformation to the system in equation (44): 


T = 


‘4 0 0 - 

0 Ac B c 

- 0 Cc D c . 


(56) 


Clearly, T is nonsingular if and only if A Z(s) has no transmission zeros at the origin. The 
transformed system has controller state equations 

x c = A C X C + B c y p (57) 

u = -y z = ~ ( C c x c + D c y p ) (58) 

where y p = (7 T ,d T ) T . Because transformation T is nonsingular, the transformed system is also 
Lyapunov-stable. The system will now be shown to be, in fact, g.a.s. 

Refer^ to figure 4(b) and see that the output y p tends to some steady-state value 
Vp — i'yssi&ss) * Because fC(s) has no zeros at the origin and is stable, its output y z = —u 
will also tend to some steady-state y z . Consequently, if y p ^ 0, the control input u will tend 
to a constant value u ^ 0. However, this contradicts the previously proven fact that u -> 0. 
Therefore, y p — » ► 0 and x c — > 0 because K{$) is stable. From LaSalle’s invariance theorem this 
proves that the system is asymptotically stable. Because V is radially unbounded, the system 
is g.a.s. Q.E.D. 

Because no assumptions were made with regard to the model order as well as to the knowledge 
of the parametric values, the stability is robust to modeling errors and parametric uncertainties. 

Remark 1: In theorem 4, if equation (47) holds with a positive-definite matrix Q c 
replacing L c L c , then the closed-loop system is g.a.s. In this case the observability and minimum 
phase property of condition 2 are not needed. 

Remark 2. The controller fC(s) stabilizes the complete plant; i.e., the system consisting 
of the rigid modes, the elastic modes, and the compensator state vector x c is g.a.s. The 
global asymptotic stability is guaranteed regardless of the number of modes in the model or 
parameter uncertainties. The order of the controller can be chosen to be any number >k. In 
other words, these results enable the design of a controller of essentially any desired order, which 
robustly stabilizes the plant. A procedure for designing £ is to choose Q z = diag (Q c , 0 k j where 
Qc — Qc ^ 0 to choose a stable A c and matrices B c and C c so that equations (47) and (48) 
are satisfied. Using equations (46) and defining 


Pz = 


Pzl 

P T 


Pz2 

PzZ 


15 



where P z \ is n c x n c and P z 3 is k x k, conditions 2 and 3 of theorem 4 can be expanded as 


PziA c + Ajp z 1 + P z^Cc + Cj PJ2 — —Qc 

P^Ac + PzzC c — 0 

bJ p z \ + dJ pT 2 = 0 

B T c P z2 + D T c P zi = I 


( 59 ) 


In addition, P z must be positive-definite. Because of the large number of free parame- 
ters ( A c ,B c ,C c ,D c ,Qc ), the use of equations (59) to obtain the compensator is generally not 
straightforward. Another method is to use the 5-domain equivalent of theorem 4. Theorem 5 
gives these equivalent conditions on K. in the 5-domain. 

Theorem 5: The closed-loop system given by equations (21), (57), and (58) is g.a.s. if JC(s) 
has no transmission zeros at 5 = 0, and IC(s)/s is MSPR. 

Proof: The proof can be obtained by a slight modification of the results of reference 29 to 
show that the theorem statement implies the conditions of theorem 4. 

The condition that JC(s)/s be MSPR is sometimes much easier to check than the conditions 
of theorem 1. For example, let JC(s) = diag[/Ci(s), . . . , ACfc(s)], where 


Ki(«) = 


*5 2 + aqjs -I- aoi 


(60) 


A straightforward analysis shows that /C(s)/s is MSPR if, and only if, k{, ctoi, Qh, 0oi, and 
are positive for i = 1, 2, . . . , k, and 

aii -Pii>0 (61) 

a li0Oi ~ OtOiPli > 0 (62) 

For higher order s, the conditions on the polynomial coefficients are harder to obtain. One 
systematic procedure for obtaining such conditions for higher order controllers is the application 
of Sturm’s theorem. (See ref. 27.) Symbolic manipulation codes can then be used to derive 
explicit inequalities. The controller design problem can be subsequently posed as a constrained 
op timiz ation problem which minimizes a given performance function. However, the case of 
fully populated IC(s) has no straightforward method of solution and remains an area for future 
research. 

The following results, which address the cases with static dissipative controllers when the 
actuators have first- and second-order dynamics, are an immediate consequence of theorem 5 
and are stated without proof. 

Corollary 5.1: For the static dissipative controller (eq. (22)), suppose that G p and G r are 
diagonal with positive entries denoted by subscript i, and actuators represented by the transfer 
function G^i(s) = ki/(s + a z ) are present in the ith control channel. Then the closed-loop system 
is g.a.s. if G r i > Gpi/ai (for * = 1 , 2 ,..., fc). 

Corollary 5.2: Suppose that the static dissipative controller also includes the feedback of 
the acceleration y a , that is, 

u = Gpyp G r y r Gaya 
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where G p , G r , and G a are diagonal with positive entries. Suppose that the actuator dynamics 
for the ith input channel are given by G^s) = ki/(s 2 + ms + i/,) with k{, m, and positive. 
Then the closed-loop system is a.s. if 

, G r i 

(* = 1 , 2 ,...,*) 

< - J ai pi 


5.3.3. Realization of K as strictly proper controller. The controller fC(s) (eqs. (57) 
and (58)) is not strictly proper because of the direct transmission term D c . From a practical 
viewpoint, a strictly proper controller is sometimes desirable because it attenuates sensor noise as 
well as high-frequency disturbances. Furthermore, the most common types of controllers, which 
include the linear quadratic Gaussian (LQG) controllers as well as the observer-pole placement 
controllers, are strictly proper; they have a first-order rolloff. In addition, the realization in 
equations (57) and (58) does not utilize the rate measurement y r . The following result states 
that fC can be realized as a strictly proper controller wherein both y p and y r are utilized. 

Theorem 6: The nonlinear plant with y p and y T as outputs is stabilized by the 

nc-dimensional controller fC given by 


x c A c%c [ B c A C L L J 


yp 

Vr 


(63) 


Vc — C c x c 

where C c is assumed to be of full rank, and an n c x k {n c > k) matrix L is a solution of 


(64) 


Da C c L — 0 (65) 

Proof: Consider the controller realization equations (57) and (58). Let 


x c = x c + Ly p (66) 

where L is an n c x k matrix. The differentiation of equation (66), use of equations (57) and (58), 
and replacement of y p with y r results in equation (63) and leads to 


y c - C c x c +(D C - C c L)y p (67) 

If L is chosen to satisfy equation (65), the strictly proper controller is given by equations (63) 
and (64). Equation (65) represents k 2 equations in kn c unknowns. If k < n c (i.e., the 
compensator order is greater than the number of plant inputs) and C c is of full rank, many 
possible solutions exist for L. The solution which minimizes the Frobenius norm of L is 

B — Cj (c c Cj'j~ 1 D c ( 68 ) 

If k = n c , equation (68) gives the unique solution L = C^ l D c ~ 

6. Linear Single-Body Spacecraft 

In this section, the case of a single-body flexible spacecraft is considered, which is a special 
case of the systems discussed in section 5. The motivation for investigating this case separately 
is that a number of spacecraft belong to this class (e.g., flexible space antennas). In addition, the 
mathematical models for this class of spacecraft are linear, which permits the use of a variety of 
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controller synthesis techniques that are not available for nonlinear systems-. However, the basic 
stability results for the linear case under static and dynamic dissipative compensation can be 
obtained as special cases of the stability results for nonlinear models. (Refer to section 5.) 


6.1. Linearized Mathematical Model 

The linearized mathematical model of the rotational dynamics of single-body flexible space 
structure can be obtained by linearizing equation (2) about the zero steady state and is given 
by _ 


where 


Mp + Kp = B 1 u 

p = (v t ,q t ) t v=(4>,0,ip) 1 


(69) 


(70) 


Here, rj represents the 3x1 attitude vector, q is the n ? xl vector of elastic modal amplitudes 
(n q = n — 3), M is the positive-definite symmetric mass-inertia matrix (note that M is constant 
in this case), K is the positive-semidefinite stiffness matrix related to the flexible degrees of 
freedom, u is the 3 x 1 input torque vector, and B = [/3,03 X n,]- The system in equation (69) 
can be transformed into modal form by using transformation Q = $p, where $ is the eigenvector 
matrix such that $ r M$ = I and 4> T K$ = C, a diagonal matrix. The resulting model is 


C + CC = $ T B T u = 


r *n n 

L*«J 


u 


(71) 


rp rjy ~ 

where $n and $ 12 are 3x3 and n q x 3 matrices, respectively, and matrix C is given by 


C = diag(0 3 , A) 


(72) 


where . 

A = diag^,W2,---,^n g ) ( 73 ) 

in which uj{(i = 1,2,..., n q ) represent the elastic mode frequencies. The first three components 
of C correspond to rigid-body modes. The rigid-body modes are controllable if, and only if, 
$ii is nonsingular. Because one torque actuator is used for each of the X-, Y-, and Z-axes, 
$n is nonsingular. Note that the mass-inertia matrix in this modal form is the 3x3 identity 
matrix. However, the model is customarily used in a slightly modified form wherein the elastic 
motion is superimposed on the rigid-body motion. (See ref. 4.) This can be achieved as follows. 
Suppose J = J T > 0 is the moment-of-inertia matrix of the spacecraft. Then, with the use of 
the transformation, £ = where the transformation matrix A is given by 


A = 


$llJ 03xn, 

. 0n,jX3 InqXriq . 


Equation (71) is transformed to 


A£ + C A£ = 

Premultiplication of the above equation by A -1 yields 


r $fl 1 

L $12 J 


U 


e + A- 1 CA^ = 


J- 1 ' 

-f-T 

$12 


u 


(74) 


(75) 


(76) 
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Note that A X CA - C and with premultiplication by diag( J, /„• ), equation (76) can be 
rewritten as 

A'z + Cz — T T u (77) 

where A — diag(J,7) and T = [/3, $12]. The inherent structural damping term B can now be 
added to give the design model as follows: 


M + B'i + ci = r T « 

(78) 

where 


A = diag(J. I nq ) 

(79) 

B = diag(0 3 , D) 

(80) 

C = diag(0 3 , A) 

(81) 

A = diag(wf,u;|,...,a;2j 

(82) 

(eqs. (81) and (82) are the same as eqs. (72) and (73), respectively), 
matrix representing the inherent damping in the elastic modes, u>j 
the elastic mode frequencies, and 

D = D > 0 is the n q x n q 
(i — 1, 2 , . . . , n q ) represents 

r=[/ 3 x 3 $ 12 ] 

(83) 

The attitude and attitude rate sensor outputs are given by 


•3 s 

11 

H 

(84) 


Vr = rz (85) 

where y p and y r are 3x1 measured position and rate vectors, respectively. 

All of the stability results of section 5 for static and dynamic dissipative controllers are 
directly applicable to this case. From theorem 2, the constant gain dissipative control law is 
given by 

u = -G p y p - G r y r (86) 

where G p and G r , which are symmetric, positive-definite, proportional and rate gain matrices, 
respectively, make the closed-loop system asymptotically stable. Furthermore from theorem 3, 
the stability is robust in the presence of monotonically nondecreasing actuator nonlinearities 
and sensor nonlinearities belonging to the (0, 00) sector. 

For this case, note that the transfer function is given by G(s) = G'(s)/s (fig. 5), where G'(s) is 
given by 


g'M-^T + E 


S 2 + 2 piLJiS + u 2 


where J is the moment-of-inertia matrix, and p t , and u t denote the rotational mode shape 
matrix, damping ratio, and natural frequency, respectively, of the ith structural mode. The 
operator represented by the transfer function G'(s) is passive. Because this operator is linear 
and time-invariant, its passivity implies that the transfer function G'(s ) is positive-real. However, 
the transfer function G(s), from u to y p , is not positive-real. This system is robustly stabilized 
by dynamic dissipative compensator £(s) satisfying the conditions of theorems 4 and 5. The 
compensator can be realized as a strictly proper one as shown in theorem 6. 
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6.2. Optimal Dynamic Dissipative Compensator 

The results of theorems 5 and 6 can be applied to check if a given model-based controller, 
such as an LQG controller, is dissipative (i.e., if it robustly stabilizes G(s)). In particular, the 
following result is obtained. 

Theorem 7: Consider the n fc th-order LQG controller given by 


q = A^q + H 


Up 

Vr 


( 88 ) 


u = Fq 

where A q is the closed-loop LQG compensator matrix 

Ac; = Aq — BqF + HC 0 


(89) 

(90) 


A$, Bq, and Cq denote the design model matrices, and F and H are the regulator and 
estimator gain matrices, respectively. This controller robustly stabilizes the system if the rational 
matrix M{s)/s is MSPR where 

M(s)= F(sl - A t ) -1 (tf i + A,H 2 )+ FH 2 (91) 

and Hi and H 2 denote the matrices consisting of the first three columns and the last three 
columns of H, respectively. 

The theorem can be proved by using the transformation x — q — H 2 y p in equation (88). 
Although a given LQG controller will not be likely to satisfy the condition of theorem 7, the 
condition can be incorporated as a constraint in the design process. The problem can be posed as 
one of minimizing a given LQG performance function with the constraint that M(s)/s is MSPR. 
Also note that theorem 7 is not limited to an LQG controller but is valid for any observer-based 
controller with control gain F and observer gain H. Another way of posing the design problem 
is to obtain the dissipative compensator which is closest to a given LQG design. The distance 
between compensators can be defined as either 

1. The distance between the compensator transfer functions in terms of H 2 or Hoo norm of 
the difference or 

2. The distance between the matrices used in the realization in terms of a matrix (i.e., 
spectral or Frobenius) norm 

For example, the dissipative compensator A c and C c matrices can be taken to be the LQG A ? 
and F matrices, respectively, and the compensator B c and D c matrices can be chosen to 
minimis 6 as follows while still satisfying the MSPR constraint: 


6 = 


B c -(Hi+A,H 2 ) 


Da - FH 2 


2 or oo 


Thus the design method usually ends up as a constrained optimization problem. 


(92) 


7. Numerical Examples 

Two numerical examples are given to demonstrate some of the results obtained in sections 4 
and 6. The first sam ple consists of a conceptual nonlinear model of a spacecraft with two flexible 
articulated appendages. The stability results for the nonlinear dissipative control law given in 
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^•^SSSSSSHSSSS 

7.1. Two-Link Flexible Space Robot 

aensity ot 2.568 x 10 kg/m^ and modulus of elasticity E = 6 34 x 10 9 b/ m 2 tu . . , 
body mass is 4030 kg and each link mass is 0 fins l™ on, • , 10 kg ' ‘ The centr al- 

central-body about loca, X, V, and “s“ y 

^^^-m^The^entral^btnlv^h 3 ^}. “ “ «“«* *•>-“« 2&£t 

Zd TSTJ^fTr ° n f bet ^ en l^^^d^^her ^ 

(quaternions) and““^ ^voTut“J^^. ** “* a,ti *“ de 

se Jnd £k « ZeTd d t d a towf b:l eam ,r th ^7?" ‘““T <™«ons; the 
purpose of si^ation, the hrst four bent^dL P Z‘ ZZ £ 

of freedom 'fo r^lf™ ngid rota,ional <*<=«««> of freedom and eight flexible degrees 

, ur for each link). The modal data were obtained from MSC/NASTRAN Thp 

mode shapes and the frequencies for links 1 and 2 are shown in figured 7 and 8 resoe^elv A 

omplete nonlinear simulation was obtained with DADS 2 , a commercially available software^ 

A rest-to-rest maneuver was considered to demonstrate the control law. The initial confitm- 

0 5 1 radTr^ar UJV f?h t0 ^{ 4 ^ rad r0tatl0n of tlie entire spacecraft about the global X-axis and 
staielf the s°7 the A reVOl y. te joint 2 ‘ Tbe <*j«*ive of the control law was to resloVe 7e zTro 

LT Because" £ ***?“ (e * (13)) Was used to -compth th" 

sdertinn 7 7, n ° kn ° Wn techni <l ues *> date for the synthesis of such controllers the 
selection of controller gams was based on trial and error. Based on several trials the fnli^’ • 

gams were found to give the desirable response: G pl = diagSo 5M lool C 7 SToT 
and G r = diaffl500 275 97 n inn ion\ * X. pl , lagtouu, 5UU, 5U0j, & p2 = diag(50, 50), 

. _ r ^ agl ° UU ’ 275 ’ 270 » 100 ’ 10 °) • As t h e system begins motion, all members move relative 
]• 7 n ° 6r ’ an ynanuc interaction exists between members. Complete nonlinear and cou 

CLf ZTo ““TheTi 1 ” The P-ame,er P r e S pou S esT s rotu°n 

fiSre 11 Th ’7 h i J ang 6 dlSplacements for the revolute joints 1 and 2 are shown in 
tiTlisnl ’ The J0mtdlsplacements deca y asymptotically and are nearly zero within 15 sec The 

L 7e f t0 gl ° bal F *’ and ap e shown in figures 12 and 3 N^e 

that the manipulator tip reaches its desired x position in about 15 sec, wheels the ^desired t, and 
^positions are reached in abo ut 35 sec. These responses 

2 I rademark of The MacNe al-Schwendler Corp., Los Angeles, CA 90041 
Trademark of Computer Aided Design Software, Inc., Oakdale, IA 52319. 
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results of section 4. 
of nonlinearities in 


The time histories of control torques are given in figures 14-17. The effects 
the model can be seen in the responses as well as in the torque profiles. 


7.2. Application to Hoop-Column Antenna 

The open-loop damping ratio is assumed to be 1 percent. A linear quadratic Gau^ian (LQG) 
controller, baLd on a design model consisting of the three rotational rigid modes and the first 
three elastic modes, was first designed to minimize 

J = lim }=£ [ (ypQpVp + Vr QrVr + U T Ruj dt (93) 

T-» 00 I J 0 v 

with O - 4 X 10 8 h Qr = 10 8 / 3 , and R = diag(0.1, 0.1,1). The actuator noise covariance 
^ensSrwl O I/3 ft-lb and the attitude and rate sensor noise covariance intensity was 
10-iCaefO 25 0 25 2.5) rad/sec and 10- 10 diag(0.25, 0.25, 2.5) rad 2 /sec 2 respectively. The 
optimal value of J was 0.6036, and the closed-loop eigenvalues for the design model an 
12th-order controller are given in table I. 

A dynamic dissipative controller, which consisted of three second-order blocks ■»_ h ‘ 
tion (60), was designed next. By using the transformation of theorem 6 with L- M) 
for JCj(s), each K,(s) can be realized as a stnctly proper controher 


Xki = 


0 1 
— ooi — Qq* 


Xki + 


Si 7 i 

foili + PliSt Si 


Vpi + Wpi 
Vri + w ri 


U =(ui,U 2,U$) T Ui — (0Qi, £*0j) x fci 


(94) 


(95) 


The constraints to be satisfied are equation (61), equation (62), and that .aw, a i»> $>*> and 
a be D ositive (i = 1 2 3). Thus, this sixth-order compensator has 18 design varia es. 

model, whereas the DDC yields guaranteed stability in the presence of higher modes as we 
parametric uncertainties. 

8. Concluding Remarks 

Stabilization of a class of nonlinear multibody flexible space systems was considf ered ’ using 
a Control laws. Robust global asymptotic stabihty can be obtained wuh 

a class ot as p centra l_ bo dy quaternion angles, relative body angles, and angular 

“££ cie wherein the centre-body motion is in the linear range 
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while all the appendages undergo unlimited motion, global asymptotic stability under a linear 
dissipative control law was proved. Furthermore, the robust stability was showed to be preserved 

in the presence of a broad class of actuator and sensor nonlinearities and a class of actuator 
dynamics. 

A class of dynamic dissipative controllers was introduced and was proved to provide global 
asymptotic stability for the case where the central-body motion is small. Dynamic dissipative 
controllers offer more design freedom than static dissipative controllers and, therefore, can 
achieve better performance and noise attenuation. 

Linear single-body spacecraft represent a special case of nonlinear multibody spacecraft, and 
therefore, the robust stability results are directly applicable to this case. 

All the stability results presented are valid in spite of unmodeled modes and parametric 
uncertainties; i.e, the stability is robust to model errors. The results have a significant practical 
value because the mathematical models of such systems usually have substantial inaccuracies, 
and the actuation and sensing devices have nonlinearities. 

Design of dissipative controllers to obtain optimal performance is, as yet, an unsolved 
problem, especially for the nonlinear case. Future work should address the development of 

systematic methods for the synthesis of both nonlinear and linear, static and dynamic dissipative 
controllers. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
December 23, 1994 
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Appendix 

Mathematical Model Derivation and System Properties 

The steps involved in the derivation of equations of motion (eq. (2)) are outlined in this 
appendix. The details of the derivation of the mass-inertia matrix M(p) and the stiffness 
matrix K are not included but can be found in reference 14. 

Kinetic Energy 

The kinetic energy of the system represented in figure 1 is given by 

T=\p T M{p)p (Al) 

where M (p) is the configuration-dependent symmetric and positive-definite mass-inertia matrix 
of the system and p is the vector of the generalized coordinates. 

Potential Energy 

Generally, the potential energy of a system has many sources (e.g., elastic displacement 
or thermal deformation). The deformations due to thermal effects are not considered in the 
formulation in this paper; however, they can easily be included in the formulation if desired. 
Thus, the potential energy is assumed to consist of the contribution only from the strain energy 
due to elastic deformation. Also, the materials under consideration are assumed to be isotropic 
in nature and to obey Hooke’s law. 

If K is the stiffness matrix of the system, then the potential energy of the system is given by 

V = l -q T Kq 

where q is the vector of the flexible degrees of freedom. If the matrix K is defined as 

°kx(n-k) 

^(n-’k)x(n-k) - 

then the potential energy of the system can be rewritten in terms of the generalized variable p 
as 

V = l -p T Kp (A2) 


K = 


Ojfcxfc 
L^n— k)xk 


Equations of Motion 

From equations (Al) and (A2), the Lagrangian of the system is formed as 

L = T-V 

For convenience, L can be rewritten in the indicial notation as 

L = T — V = ^Mijpipj - V (A3) 

ij 

The Euler-Lagrange equations for the system can then be derived from 



24 



where F k axe generalized forces from the nonconservative force field. From evaluation of the 
derivatives, 

§^=Y, M kjPj (A5) 

3 


and 


i 3 


^kjPj + y~! 


9M kj . . 




dp 


— PiPj 


Also 


dL _ 1 ^ 


dpk 2 dp k PlPj dp k 

l yJ 

Thus, the Euler-Lagrange equations can be written as 

w .. , v 1 dM x j \ 5V" 


dV 


(k = 1,2, ...,n) 


By interchanging the order of summation and taking advantage of symmetry, 


(A6) 


(A7) 


(A8) 


( dM k j \ . . 9M kl dM ki \ . . 

£( * = - 2 Z{ + ^')r* 


(A9) 

Hence, 

T (dM kj 1 dMij\ ^1 (dM kj d M kt dM {j 

ij ' 2 ®Pk J 1 3 2 \ dpi dpj dp k 

^ PiPj 

(A10) 

The terms 

c m - l -( aM « + aM * 

2 \ dpi dpj dp k ) 


(All) 

are known as Christoffel symbols. For each fixed k, note that = Cj ki . 

Also, 


dV ^ 
dp k ~ KkjPj 


(A12) 

Finally, the Euler-Lagrange equations of motion can be written as 



Y1 M kjPj + C ijkPiPj + D kjP + K kjPj = F k (k= 1,2, 

...,n) 

(A13) 


3 »J 

where D is the inherent structural damping matrix and Dp is the vector of nonconservative 
forces. 

Of the four terms on the left side of equation (A13), the first term involves second derivatives 
of the generalized coordinates p. The second term consists of centrifugal terms (e.g., pf) and 
Coriolis terms (e.g., pipj, i ± j). In general, the coefficients C tJ are functions of p. The third 
term involves only the first derivatives of p and corresponds to the dissipative forces due to 
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inherent damping. The fourth term, which involves only p, arises from the differentiation of the 
potential energy. 

In matrix- vector notation, equation (A13) can be written as 

M(p)p + C{p,p)p + Dp + Kp = F (A14) 


The k,jth element of the matrix C(p,p) is defined as 

n 

Ckj = ^ ] Cijk(p)fii 


i=l 


1 { 9M k j dMkj — dM t j \ . 

~ 2^2\ dpi dpj dp k ) Vl 


(A15) 


An important property of the systems whose equations of motion are given by equation (A14) 
is derived next. This property is pivotal to the stability results obtained in sections 4 and 5. 

Theorem Al: The matrix S = M (p ) - 2C(p,p)is skew-symmetric. 

Proof: The k, jth element of the time derivative of the mass-inertia matrix M(p) is given 
by the chain rule as 


^ kj — yi 


i=i 


dpi 


-Pi 


(A16) 


Therefore, the k,jth component of S = M — 2 C is given by 

$kj — Mkj — 2 Ckj 


= E 


\9M kj 

j a y + 

. dpi 

V dpi 

( dM{j 

9M k A. 

V dpk 

dpj / 


dp j 


dpk 


')]* 


(A17) 


Because the inertia matrix is symmetric (i.e., Mjj = Mji), the interchange of the indices k 
and j in equation (A17) results in 

Sjk = ~ S kj ( A18 ) 

This completes the proof. Q.E.D. 

Theorem Al can be used to prove that the system given by equations (2) has the important 
property of passivity as defined in reference 3. 

Theorem A2: The input-output map from u to y T is passive, i.e., with zero initial conditions, 


[ T yj(t)u(t)>0 (VT > 0) 

J 0 


(A19) 


for all u(t) belonging to the extended Lebesgue space L Re- 
proof: Premultiplication of both sides of equations (2) by p T and integration result in 

J ^p T M(p)p + p T C{j),p)p + p T Dp + p T Kp^ = yjudt (A20) 
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Note that 


J t [p T M (p)p] = 2 p T M(p)p + p T M(p)p 
Application of theorem Al and simplification yields 

^P T (T)M\p(TIP(T) + £ pT D p dt + \p T {T)Kp{T) = £ y T u dt 

Because the left side of equation (A22) is nonnegative for ail T > 0, this gives the required result. 

Q.E.D. 


(A21) 


(A22) 
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Table I. Closed-Loop Eigenvalues for LQG Controller 


Regulator 

Estimator 

-0.0238 ± 0.05421 
-0.0721 ± 0.0964i 
-0.0754 ± O.lOOOz 
-0.3169 ±0.8108i 
-0.2388 ± 1.3554i 
-0.3375 ± 1.7030z 

-0.0240 ± 0.0544Z 
-0.0720 ± 0.0965i 
-0.0725 ± 0.0963z 
-0.4050 ± 0.7703 z 
-0.3334 ± 1.3333z 
-0.5104 ± 1.6562z 


Table II. Closed-Loop Eigenvalues for 
Dynamic Dissipative Controller 

-0.0035 ± O.Ulvm 
-0.0183 ± 0.0458z 
-0.0160 ± 0.0502z 
-0.3419 ± 0.5913z 
-0.7179 ± 0.6428z 
-0.8479 ± 0.5653z 
-0.6482 ± 1.6451z 
-0.4536 ± 2.1473i 
-0.3764 ± 2.5522z 
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Branch 



Figure 1. Multibody system 



(a) (0, oc)-sector nonlinearity. 



(b) (0, oo)-sector monotonically nondecreasing nonlinearity. 
Figure 2. Examples of nonlinearities. 
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Figure 3. Feedback closed-loop configuration. 



(a) Compensator with rate input. 



(b) Compensator with position input. 
Figure 4. Rearrangement of feedback loops. 
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Figure 6. Flexible space robot. 
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Displacement Displacement 



(a) Mode 1; frequency = 2.14 Hz. 


(b) Mode 2; frequency = 8.58 Hz. 



(c) Mode 3; frequency = 19.3 Hz. (d) Mode 4; frequency = 34.32 Hz. 

Figure 7. Mode shapes of link 1. 
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Displacement 



0 5 10 15 20 0 5 10 15 20 


Node number Node number 

(a) Mode 1; frequency = 3.3 Hz. (b) Mode 2; frequency = 10.8 Hz. 



(c) Mode 3; frequency = 22.5 Hz. (d) Mode 4; frequency = 38.3 Hz. 

Figure 8. Mode shapes of link 2. 
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Figure 9. Euler parameter qj. 


1.0 x 10 -3 



Figure 10. Euler parameters a 2 and a 3 
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Figure 11. Revolute joint angular displacement. 
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Torque, N-m 




Figure 15. Control torques «2 and «3- 
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Figure 17. Control torque u$ at revolute joint 2 
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First 10 natural frequencies, rad/sec: 
0.75, 1.35, 1.70, 3.18,4.53, 
5.59, 5.78, 6.84, 7.40, 8.78 



Figure 18. Hoop-column antenna concept. 
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